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ABSTRACT 

Using resistive magnetohydrodynamics simulation, we investigate circumstel- 
lar disk formation in a strongly magnetized cloud. As the initial state, an isolated 
cloud core embedded in a low-density interstellar medium with a uniform mag- 
netic field is adopted. The cloud evolution is calculated until almost all gas inside 
the initial cloud falls onto either the circumstellar disk or a protostar, and a part 
of the gas is ejected into the interstellar medium by the protostellar outflow 
driven by the circumstellar disk. In the early main accretion phase, the disk size 
is limited to ~10 AU because the angular momentum of the circumstellar disk is 
effectively transferred by both magnetic braking and the protostellar outflow. In 
the later main accretion phase, however, the circumstellar disk grows rapidly and 
exceeds ^ 100 AU by the end of the main accretion phase. This rapid growth of 
the circumstellar disk is caused by the depletion of the infalling envelope, while 
magnetic braking is effective when the infalling envelope is more massive than 
the circumstellar disk. The infalling envelope cannot brake the circumstellar disk 
when the latter is more massive than the former. In addition, the protostellar 
outflow weakens and disappears in the later main accretion phase, because the 
outflow is powered by gas accretion onto the circumstellar disk. Although the 
circumstellar disk formed in a magnetized cloud is considerably smaller than that 
in an unmagnetized cloud, a circumstellar disk exceeding 100 AU can form even 
in a strongly magnetized cloud. 

Subject headings: accretion, accretion disks: ISM: clouds — stars: formation — 
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Introduction 



The magnetic field is an essential ingredient in present-day star formation. Observations 
have shown that mol ecular cloud cores have magnetic energy compara ble to gravitational en- 
ergy JCrutcherlll999h. The cloud cores also have angular momentum (lArquilla Sz Goldsmith 



1986 



Goodman et al. I Il993l ; ICaselli et al.l 120021 ) . Stars are born in such cloud cores. As 



the cloud collapses, the central region of the collapsing cloud rotates rapidly, and this rapid 
rotation suppresses further collapse and subsequent protostar formation. In the collapsing 
cloud, however, the magnetic field can effectively transfer the excess angular momentum from 
the center of the cloud by magnetic braking (and protostellar outflow), and promote further 
collapse and protostar formation. Moreover, after the protostar formation, the magnetic 
field greatly contributes to the evolution of a protostar and a circumstellar disk. Because 
the angular momentum in the circumstellar disk is transferred by magnetic effects, the gas 
in the circumstellar disk can fall onto the protostar, increasing the protostellar mass in the 
main accretion phase. Thus, it may be expected that only a small-size (or no) disk ap- 
pears around the protostar if the angular momentum of the disk is effectively transferred 
by magnetic effects. In various star-forming regions, however, many circumstellar disks wit h 
a size of several hundred AU can be observed (e.g., ICalvet et al.l I200CH ; iNatta et al.l 120001 ). 
indicating that some of the angular momentum survives the magnetic braking catastrophe 
and contributes to the growth of the circumstellar disk. 

The formation and evolution of the circumstellar disk are important not only for star 
formation but also for planet formation. The size and mass of the circumstellar disk deter- 
mi ne the mode of forma tion of a planet (i.e., core accretion or gravitational instability mode; 
see iDurisen et al.l 120071 ). The angular velocity of the molecular cloud core that is the host 



cloud for star formation creates a circumstellar disk in the star formation process. Thus, to 
understand the formation of the circumstellar disk, we need to investigate the evolution of 
the molecular cloud core from the prestellar stage. Numerical simulations can reproduce the 
disk formation in the molecular cloud core. However, because there is a large difference in 
both spatial and time scales between the molecular cloud core and the circumstellar disk, 
we need special numerical techniques such as adaptive mesh refinement or nested grid and 
smoothed particle hydrodynamics. Using these techniques, several studies have investig ated 
the disk forma t ion in the molecular cloud core from the prestellar stage. iBatd ( 120101 ) and 
Machida et al.l ( 12 01 Obi ) investigated the formation of the protostar and circumstellar disk in 
an unmagnetized cloud; they showed that the first core that forms before protostar formation 
in the g as-collapsing phase direc tly evolves into a circu mstellar disk after the protostar for- 



mation. 



Walch et al.l ( 120091 ) and iMachida et al.l ( 12010al ) showed that the circumstellar disk 



which originates from the first core, is more massive than or comparable to the protostar 
in the main accretion phase and has a size of several hundred AU when the host cloud is 
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unmagnetized. 



Using numerical simulations, the disk formation in magnetized clouds has also been 
investigated. These studies reported that the calculated size of the circumstellar disk ap- 
pearing in a magnetized c loud i s much smaller than the observed sizes. Using an ideal MHD 
calculation, iMellon fc Lil ( 120081 ) showed that the circumstellar disk forme d in a magnetized 
cloud has a size of ~ 10 AU at maximum. iHennebelle fc Ciardil (120091 ) also showed that 
large-size disk formation is suppressed in a magnetized cloud. This is because the angular 
momentum of the circumstellar disk (or infalling gas onto the circumstellar disk) is exces- 
sively transferred outward by magnetic braking, and the gas in the disk effectively falls onto 
the protostar. In reality, however, the magnetic field dissipates in a high-density gas region 
by Ohmic dissip ation and ambipolar diffusion, because the degree of ionization is consid- 
erably low (e.g., iNakano et al. 1 120021 ) . Such a high-density gas region, where the neutral 



gas is decoupled from the magnetic field, corresponds well to the disk- forming region in th e 
early main accretion phase (~ 0.1 — 10 AU, iMachida et al.ll2007l l2008at llnutsuka et al.ll2010l ) . 
Thus, if the magnetic field can sufficiently dissipate around the protostar, magnetic brak- 
ing becomes ineffective and a large-size circumstellar disk comparable to observations may 
appear even i n stro ngly magnetized cloud s . In c ontrast, with non ideal MHD calculations, 
Mellon &: Lil ( 120091 ) and iDuffin & Pudritzl (120091 ) reported that the disk size is considerably 
smaller than observations even when the magnetic field dissipates in the high-density gas 
region, although the disk in non ideal MHD calculations is somewhat larger than that in 
ideal MHD calculations. 

These studies, however, overlooked the effect of the finite mass of the host cloud for 
star and circumstellar disk formation. In the early main accretion phase, a massive gas 
envelope (or remnant of the molecular cloud) remains around the circumstellar disk. The 
massive envelope rotates much slower than the circumstellar disk, and thus it can brake 
the circumstellar disk through the magnetic field lines that connect the less massive rapidly 
rotating disk to the more massive slowly rotating envelope. In other words, the angular 
momentum of the circumstellar disk is transferred into the infalling envelope by torsional 
Alfven waves. This is the magnetic braking mechanism for the circumstellar disk formation. 
Thus, even with very strong magnetic fields, the circumstellar disk rotation velocity does not 
decrease when no massive envelope exists around the disk. The mass of the infalling envelope 
decreases with time in the main accretion phase when an isolated molecular cloud core is 
assumed. The infalling envelope disappears by the end of the main accr etion phase: a part 
of th e envelope is blown aw ay from the cloud by the protostellar outflow (IMatzner &: McKee 
2000l ; IMachida et al.l 120091 ) , and the remainder falls onto either the circumstellar disk or 
the protostar. Thus, we can expect that magnetic braking becomes gradually ineffective 
as the mass of the infalling envelope decreases, and it rarely works when the mass of the 
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circumstellar disk is much larger than that of the inf ailing envelope. Note that, in reality, 
magnetic braking can work even after the infalling envelope disappears, because the low- 
density interstellar medium can receive angular momentum from the circumstellar disk. 
Also note that because the efficiency of magnetic braking decreases as the ambient density 
decreases, the angular momentum of the circumstellar disk is not effectively transferred into 
the low-density interstellar medium (see, §1]). 

In addition, in an isolated cloud core, although magnetic braking can redistribute the 
angular momentum inside the cloud, it cannot effectively transfer the angular momentum 
outside the cloud (or into the interstellar medium). The infalling envelope receives the 
angular momentum from the circumstellar disk by magnetic braking, while it finally falls 
onto the circumstellar disk with the angular momentum by the end of the main accretion 
phase. However, when the infalling envelope is depleted and magnetic braking becomes 
ineffective, no sufficiently large disk appears if insufficient mass to form the disk remains in 
the infalling envelope. 

In summary, to investigate the real size and mass of the circumstellar disk, we have to 
investigate the disk evolution in the collapsing cloud core, at least till the end of the main 
accretion phase (or until the infalling gas disappears). In previous studies, the calculation 
was stopped much before the infalling envelope disappears. Thus, it is natural that mag- 
netic braking is effective and no large-size disk appears during this phase (i.e., early main 
accretion phase). In this study, however, we calculated the evolution of the circumstellar 
disk until almost all gas of the infalling envelope disappears and found that a large-size disk 
comparable to observations appears when the infalling envelope becomes less massive than 
the circumstellar disk. The structure of this paper is as follows. The framework of our 
models and the numerical method are given in §2. The numerical results are presented in 
§3. We discuss the magnetic braking timescale and the effect of the infalling envelope on the 
circumstellar disk formation in 54. 



To study the formation and evolution of a circumstellar disk in a magnetized molecular 
cloud core, we solve the three-dimensional resistive MHD equations including self-gravity: 



2. Model Settings 



2.1. Basic Equations and Initial State 



dp 
dt 



+ V ■ (pv) = 0, 
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OB 

~dt 



V x (v x B) +rjV 2 B, 
V 2 = 4ttG P , 



(3) 
(4) 



where p, v, P, B, 77, and <fi denote density, velocity, pressure, magnetic flux density, resistiv- 
ity, and gravitational potenti al, respectively. To mimic the thermal evolution calculated by 



Vorobyov fc Basu 



2006 



Masunaga fc Inutsuka ( 2000T). we adopt the piece wise polytrop i c equ a tion of state (see also 



Machida et all 120071 : Mellon fc Lill2008[ 120091 : llnutsuka et al.ll2010h 



as 



c loP 



1 + 1^ 



where c Sj o = 190 ms 1 and p c = 3.84 x 10 14 g cm 



2/3' 



(5) 



10 10 cm- 3 ). With equation © 



the gas behaves isothermally for n <^ 10 cm" and adiabatically for n ^ 10 cm" . For a 
realistic evolution o f the magnetic f i eld in the circumstellar disk, we adopt the resistivity 77 as 
the fiducial value in lMachida et al.l ( 120071 ). in which the Ohmic dissipation be comes effective 



for 10 11 cm 3 <^ n <^ 10 15 cm 3 (for details, see Eqs. [9] and [10], and Fig. 1 of Machida et al. 



2007J). 



As the initial state, we take a spherical cloud with a Bonnor-Ebert (B E) density profile 



which is widely rec ognized as the density profile of a prestellar cloud (lAlves et al.l 12001 
Kandori et al.l 120051 ) . The initial cloud extends up to r = 2 ?"be; where t"be is the critical BE 



radius. Outside the initial cloud (r > 2 tbe), a uniform density is adopted to mimic the static 
interstellar medium. We a dd nonaxisymmetric de nsity perturbation of the m = 2 mode to the 
initial core (for details, see iMachida et al.ll2010al ). This perturbation ensures that the center 
of gravity is always located at the origin. We ignore gravity (self-gravity and gravity of the 
protostar) in the region of r > 2 the- In addition, we prohibit the gas inflow at r = 2tbe in 
order to suppress mass input from the interstellar medium into the gravitationally collapsing 
cloud core. Thus, only the gas inside r < 2 tbe collapses to form the circumstellar disk and 
protostar. Note that although the protostellar outflow driven by the circumstellar disk 
propagates into the region of r > 2tbe and disturbs the interstellar medium over time, 
we can safely calculate the mass accretion process onto the circumstellar disk because the 
outflow propagating into the interstellar medium does not affect the inner cloud region 
(r < 2 tee)- For the BE density profile, we adopt the central density as n CiQ = 3 x 10 4 cm -3 , 
3 x 10 6 cm -3 (fiducial model), and 3 x 10 7 cm -3 with an isothermal temperature of T = 10 K. 
To promote contraction, we increase the density by a factor of /=1.68, where / is the 
density enhancement factor that represents the stability of the initial cloud. With / = 1.68, 



0.5 for r < tbe and 0.42 for r < 2 tbe, where and E, 



grav 



the initial cloud has Eth/l-^grav 
are the thermal and gravitational energy, respectively. Although the central density of the 
initial cloud has values of n Cj o = 3 x 10 4 / cm -3 , 3 x 10 6 / cm -3 and 3 x 10 7 f cm -3 , we describe 



-6 - 



the initial cloud density without the density enhancement factor / as n c p = 3 x 10 4 cm -3 in 
the following. Note that because the critical BE radius (or initial cloud size) is determined 
by the central density without /, the initial density without / is a more suitable parameter 
for describing the initial cloud. The density and enclosed mass against the initial cloud 
radius are plotted in Figured] for the model with n Ci o = 3 x 10 6 cm~ 3 and 3 x 10 7 cm~ 3 . The 
density contrast between the center of the cloud (r = 0) and cloud boundary (r = 2 tbe) is 
p(r = 0)/p(r = 2t B e) = 82. The cloud with a central density of n c>0 = 3 x 10 6 cm" 3 (fiducial 
model) has a radius of r c \ = 5400 AU and mass of M c \ = 1M . On the other hand, the 
cloud with n c .o = 3 x 10 4 cm" 3 and 3 x 10 7 cm" 3 have a radius of r c \ = 2 tbe = 5.4 x 10 4 AU 
and 1740 AU, and a mass of M c \ = 10 M Q and 0.3 M Q , respectively. As seen in Figure (TJ 
because we adopted the BE density profile, the density has a nearly flat distribution around 
the center of the cloud, while it drops in proportion to r~ 2 in the outer envelope. 

The cloud rotation generates a rotating disk around the protostar as the cloud collapses. 
However, it is considered that the magnetic field suppresses the formation of the rotationally 
supported disk in the collapsing cloud, because the magnetic field effectively transfers the 
angular momentum outward from the center of the cloud. In this study, to investigate the 
disk formation in an initially rotating magnetized cloud, we parameterized the magnetic field 
strength and rotation rate of the initial cloud. We adopt a uniform magnetic field along the 
z-axis in the whole computational domain and rigid rotation along the z-axis in the region 
of r < 2tbe- Thus, the rotation axis is parallel to the magnetic field lines at the initial state. 
To specify the initial cloud, we use two parameters describing the magnetic field strength 
(a) and rotation rate (go). They are normalized by the central quantities and described as 

r2 

where B z is the magnetic flux density, and 

where f2 C) o is the angular velocity at the center of the cloud. Model names and their param- 
eters are listed in Table [TJ We adopt a = 0.1, 0.3, and 1, and u = 0.05, 0.1 and 0.15, as 
described in Table [TJ 

In this study, we investigate disk formation until almost all the mass of the initial cloud 
falls onto either the circumstellar disk or the protostar. Note that a part of the cloud mass 
is ejected from the collapsing cloud by protostellar outflow. Thus, quantities describing the 
whole cloud may be useful to specify our models. The parameter a (or the magnetic flux 
density at the center of the cloud) can be connected to the mass-to-flux ratio of the whole 
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cloud. The mass-to-flux ratio is described as 

M(r) M(r) 



$(r) 7rr 2 5 

where M(r) is the mass contained within the radius r and $(r) is the magnetic flux threading 
the cloud. There exists a critical value of M/$ below which a clo ud is supported against 



gravit y by the magnetic field. For a cloud with uniform density, IMouschovias fc Spitzer 



( 119761 ) derived a critical mass-to-flux ratio 

(9) 



$ / „,= 3tt V G 



where the constant ( = 0.48 (ITomisaka et al. lll988al Jbl). The mass-to-flux ratio normalized 



by the critical value A is described as 

The mass-to-flux ratios normalized by the critical value for different models (or different 
a) are plotted against the enclosed mass (left) and the radius (right) of the initial cloud in 
the upper panels of Figure El Reflecting the BE density profile, the mass-to-flux ratio has 
a radial distribution inside the initial cloud. For example, with a — 1, the mass-to-flux 
ratio has a peak at ~ 0.1 M (~ 0.05 M©) or ~ 10 3 AU (~ 400 AU) for the model with 
n c ,o = 3 x 10 6 cm~ 3 (n C) o = 3 x 10 7 cm~ 3 ). In addition, with a = 1 and n Cj o = 3 x 10 6 cm -3 , 
the initial cloud has a mass-to-flux ratio of A(5400AU) ~ 1 as the whole. Hereafter, for 
simplicity, we describe the mass-to- flux ratio of the whole cloud (r < 2 tbe) as A (= A[2 tbe)]), 
and call it the mass-to-flux ratio. As listed in Table [Q, we adopt A = 1, 2, 3 (magnetized 
model) and oo (unmagnetized model). The observations indicate that molecular cloud cores 
have the mass- to-flux ratio in the range of 0.8 <^ A <^ 7.2 with a median value of A ~ 2 



(ICrutcherlll999l ). Thus, the mass-to-flux ratios adopted in this paper are within a reasonable 



range. 

As the parameter for the rotation of the whole cloud, we define the initial ratio of the 
rotational to the gravitational energy as 

00 - ifr, (11) 

where E T is the rotational energy of the initial cloud in the region of r < 2 tbe- As listed in 
Table [U models have the parameter in the range of 0.007 < (3q < 0.04. Observations showed 
that molecular cloud cores have the rotational energy in the range of 10~ 4 ^ flo ^ 0.07 with a 
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median value of /3q ~ 0.02 (IGoodman et al. Ill993l ; ICaselli et al.l 120021 ). Thus, the rotational 
energies in our models are also within reasonable ranges, compared with the observations. 
In addition, we also listed the centrifugal radius of each model in Table [TJ The centrifugal 
radius is defined as 

j( r ) n 



^ = 1TWY (12) 

where j(r) is the specific angular momentum. The centrifugal radii for models with different 
u or different n Ci0 are plotted against the enclosed mass (left) and the radius (right) of the 
initial cloud in the lower panels of Figure [2j In the figure, the centrifugal radius gradually 
increases as the enclosed mass increases, because the gas far from the center of the cloud 
has a large specific angular momentum. As listed in Table [U models have the centrifugal 
radius in the range of 540AU < r c f < 2.15 x 10 4 AU as the whole cloud (r < 2tbe)- Thus, 
when the angular momentum is well conserved in the collapsing cloud, a disk with a size of 
^ 1000 AU can form around the protostar. 



sec 



2.2. Numerical Method and Boundary Condition 

For calculations on a l arge spatial scale, the nested grid method is adopted (for details, 



Machida et al. I l2005al lbl). Each level of a rectangular grid has the same number of 
cells (64 x 64 x 32). The calculation is first performed with five grid levels (I = 1 — 5). 
The box size of the coarsest grid / = 1 is chosen to be L\ = 2 5 rBE, which corresponds to 
~ 8.7xl0 5 AU (n Ci0 = 3xl0 4 cm~ 3 ), ~ 8.7xl0 4 AU (n c>0 = 3xl0 6 cm" 3 ) and ~ 2.7 x 10 4 AU 
(pr.,n = 3 x 10 7 cm~ 3 ). A new finer grid is generated before the Jeans condition is violated 



( jTruelove et al. 1119971 ). For model with n Cj o = 3 x 10 6 cm -3 , the maximum level of grids is 
l m£H = 12, which has a box size of L\2 = 43 AU and a cell width of A; max = 0.67 AU. Models 
with different cloud sizes (or different initial central densities) have almost the same spatial 
resolution (A/ max ~ 0.5 AU) of the maximum level grid, in which the maximum grid level is 
different in models with different initial central density. 

As described in §2.1^ we adopt gravity only in the region of r < 2 tbe to mimic the star- 
forming cloud core embedded in the interstellar medium, while a large domain (r > 2tbe, 
interstellar medium) is left behind from the gravitational contraction. Thus, in our setting, 
the star forming cloud with a size of 5700 AU (5.7 x 10 4 AU, 1740 AU) exists in the interstellar 
medium with a box size of ~ 8.7 x 10 4 AU (~ 8.7 x 10 5 AU, ~ 2.7 x 10 4 AU) for a model with 
n c ,o = 3 x 10 6 cm -3 (3 x 10 4 cm -3 , 3 x 10 7 cm -3 ). The purpose of this study is to investigate 
the effect of magnetic braking on disk formation. The angular momentum is transferred by 
magnetic braking with the Alfven speed. To suppress artificial reflection of Alfven waves 
at the computational boundary, we adopt a large simulation box to reproduce a realistic 
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interstellar medium around the star forming cloud. The Alfven speed in the interstellar 
medium (r > 2 tbe) for the model with the strongest magnetic field (B zQ = 176 fiG) and 
n c = 3 x 10 6 cm -3 is va = 1.3 kms -1 . Thus, for the Alfven wave to reach the computational 
boundary from the center of the cloud, it takes 3 x 10 5 yr (= Li/ua), which corresponds 
to about 30 times the free-fall timescale (tg iC = 1.0 x 10 4 yr), where the freefall timescale 
tfj jC is normalized by the central density in the initial cloud. In all models, we stopped the 
calculation within t < 20tg jC . Thus, the Alfven wave generated at the center of the cloud 
(or the computational boundary) never reaches the computational boundary (or the center 
of the cloud). With this setting, we can correctly treat magnetic braking in the collapsing 
cloud. 

In the collapsing cloud core, we assume the protostar formation to occur when the 
number density exceeds n > 10 12 cm~ 3 at the center of the cloud. To model the protostar, 
we adopt a sink around the center of the computational domain. In the region of r < 
'"sink = 1 AU, gas having a number density n > 10 12 cm -3 is removed from the computational 



domain and added to the protostar as gravity in each timestep (for details, see iMachida et al. 



20091 . l2010af ). This treatment of the sink makes it possible to calculate the evolution of the 
collapsing cloud and circumstellar disk for a longer duration. In addition, inside the sink, 
the magnetic flux is removed by Ohmic dissipation, because such a region has a magnetic 



Reynolds Re number exceeding unity Re > 1 (for details, see IMachida et al.l 120071 ). 



3. Results 

To understand the process of formation of a circumstellar disk in a magnetized cloud, 
we prepared nine different initial clouds with various masses, radii, magnetic field strengths, 
and angular velocities as listed in Table [TJ These clouds have masses of M c i i0 = 0.3 — 10 M Q 
and radii of R c i,o = 1740 — 5.4 x 10 4 AU. The magnetic field strengths Bq and angular 
velocities VLq at the center of the initial cloud are in the range of B = 17.6 — 566 fiG and 
Qq = (2.3 — 74) x 10~ 14 s -1 , respectively. In the following subsections, after we describe the 
disk formation for a typical model (model 1), we compare the disk formation in models with 
different initial settings. 



3.1. Disk Formation for a Typical Model 



In this subsection, we describe the disk formation for model 1, which has the parameter 
of A = 1. Figure [3] shows the time sequence of the density distribution around the center of 
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the collapsing cloud on the equatorial plane for model 1. The elapsed time in units of the 
initial free-fall timescale tg tC and year after the calculation begins are described in each panel. 
The protostellar, circumstellar disk, and envelope masses (M ps , M disk , and M env , respectively) 
are also described in each panel. The protostellar mass M ps is defined as the mass falling 
into the sink. To determine the circumstellar disk in the infalling envelope, we estimated 
the velocity ratio R v (= \v r /v^\, where v r and are the radial and azimuthal velocity, 
respectively) in each cell, and specified the most distant cell having R v < 1 from the center 
of the cloud. Then, we defined the disk radius r d as the distance of the cell furthest from the 
origin, and the disk boundary density pd,b as the density of the most distant cell. Finally, we 
defined the circumstellar disk as a region that h as a density of p > pd,b and estimated the 



circumstellar mass Mdi S k- iMachida et al.l (j2010al ) showed that the Keplerian rotating disk is 



well identified with this criterion. In addition, in each timestep, we estimated the total mass 
M tot that is sum of the total mass inside the gravitational sphere inside r < 2 r B E and the 
protostellar mass. Then, we subtracted the circumstellar disk mass and protostellar mass 
from the total mass and derived the mass of the infalling envelope as 

M env = Mtot - M disk - M ps . (13) 

We confirmed that the sum of the protostellar mass and the total mass is well conserved 
during the calculation. Note that, in the later main accretion phase, because the mass 
is ejected from the collapsing cloud by the protostellar outflow, mass conservation inside 
r < 2 tee is not realized. 

Figure [3] shows that the circumstellar disk gradually increases its size with time. Fig- 
ure |3]a indicates that t he circumstellar disk has a size of ~ 10 AU in the early main accretion 



phase, as described in iMellon &: Lil (120081 120091 ). At this epoch, the circumstellar disk has 
a mass M^sk = 0.07 5 Mp, which is about twice the protostellar mass (M ps = 0.038 M ). 



Inutsuka et al.l (120101 ) showed that the circumstellar disk mass is greater than the protostel- 
lar mass in the early main accretion phase, because the first adiabatic core, which is formed 
before the protostar formation with a mass of ~ 0.1 — 0.01 M Q , directly evolves into the 
circumstellar disk after the protostar formation. In addition, at the same epoch, the mass of 
the infalling envelope is much larger than that of the circumstellar disk (M env > lOM^^). 

In Figure HI the radial distribution of the azimuthal velocity normalized by the Keplerian 
velocity (a), ratio of radial to azimuthal velocity (6), surface density (c) and plasma beta (d) 
at each epoch corresponding to each panel in Figure [3] are plotted. The azimuthal velocity 
and ratio of the radial to azimuthal velocity are azimuthally averaged on the equatorial 
plane. The surface density is estimated as 

a e (r,9)= [ C " p(r,9,z)dz, (14) 
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where z cr \ = 50 AU is adopted, and azimuthally averaged. Note that the surface density 
depends only a little on z cri when z cri ^> 10 AU (and/or z CTi ^ 100 AU), because the disk 
mass is concentrated on the equatorial plane. The plasma beta is defined in each cell on the 
equatorial plane as 

8ttP , s 

P = -#r, (15) 

where B and P are the magnetic flux density and thermal pressure on each cell, respectively, 
which are also azimuthally averaged on the equatorial plane. 

At the same epoch as Figure [3]a (t = 3.94 x 10 4 yr, black lines in Fig. H]), the disk has a 
Keplerian velocity in the range of r ^ 10 AU (Fig. H]a), while the disk structure extends up 
to ~ 20 - 40 AU (Fig. 1c). In the range of 10 AU £ r £ 40 AU, the disk has 20 - 80% of the 
Keplerian velocity. Figure H]& shows that the azimuthal velocity is greater than the radial 
velocity in the range of r ^ 20 AU. Thus, the disk is mainly supported by the centrifugal 
force in the range of r ^ 20 AU, while it is also supported by the thermal pressure gradient 
and Lorentz force in addition to the centrifugal force in the range of 20 AU <; r <^ 40 AU. 
The black line in Figure shows the rapid increase of the plasma beta at r ~ 30 AU. 
Because the thermal pressure increases in the high-density region, the plasma beta tends to 
increase inside the disk. In addition, the magne tic field dissipates by Ohmic dissipation in 



such a high- density region (INakano et al. II2002I ). 



Figures [3] and H] indicate that the circumstellar disk increases its size with time. At 
t = 7.67 x 10 4 yr (Fig. [3]fr and red lines in Fig. HJ), the circumstellar disk extends up to 
~ 60 AU and has a mass of ~ 0.29 M & . Figure H]a shows that the disk rotates with a 
Keplerian velocity in the range of r ^ 50 — 60 AU at this epoch. The radial distribution of 
the surface density indicates that the disk extends up to r ~ 50 — 70 AU. In addition, the 
plasma beta rapidly increases at r ~ 60 AU and reaches /3 p ~ 10 6 near the origin. Thus, the 
whole region of the disk (r <^ 60 AU) is mainly supported by the centrifugal force. In other 
words, the Keplerian rotating disk extends up to ~ 60 AU at this epoch. 

The circumstellar disk mass becomes comparable to the mass of the infalling envelope 
at t ~ 10 5 yr, which corresponds to t ~ 10% jC (Fig. [3]c and blue lines in Fig. H]). At the 
initial state, the density at the cloud boundary (i.e., r = 2 r be) is 82 times lower than that 
at the center of the cloud (r = 0). Thus, the freefall timescale at the cloud boundary tg^ is 
about nine times longer than that at center of the cloud £ff )C (i.e., tg^ = y/82fs jC = 9.1 fg c ). 
Thus, we can expect that a large fraction of the initial cloud mass has already fallen onto the 
central region at t 3> 9/ff, c - Note that the cloud r otation and magnetic field sli ghtly slows 



the cloud collapse and accretion onto the disk (e.g.. iMachida et al. Il2005all2006al ). FigureEJc 
shows that the size of the circumstellar disk exceeds 100 AU and the spiral structure appears 
in the disk at t = 1.06 x 10 5 yr = 10.2 tff c . The blue lines in Figure H] also indicate that the 
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Keplerian rotating disk extends up to r ^ 100 AU at this epoch. 



At t = 1.25 x 10 5 yr = 12.04 % c (Fig. [3H and green lines in Fig. HJ), the circumstellar 
disk mass (M^k — 0.49 M ) is much greater than the mass of the infalling envelope (M env = 
0.23 M Q ). Thus, as described in ^TJ it is expected that magnetic braking is not effective at 
this epoch (for details, see §4]). The green lines in Figure H] indicates that the circumstellar 
disk with a Keplerian velocity has a size of ~ 200 AU (see Fig. |Uc) at this epoch. Thus, the 
circumstellar disk exceeding 100 AU can form even in a strongly magnetized cloud core. 

Figure [5] shows the time sequence of the density distribution for model 1 at the same 
epochs as in Figure [3] but for the y = plane. The outflow is driven near the cir- 
cumstellar disk, as seen in Figure [5]a-c. This type of outflow was seen in many previ- 



ous studies fe.g.. iTomisakal 120021: Machida et al. 112004 12 005b; Machida et al.ll2006l . 12008b 



Hennebelle fc Fromang!l2008l ; IPuffin k, Pudritj 2009 ; Tomida et al.lboiol ). In each panel, the 



Keplerian rotation disk is indicated by the black circle, inside which the gas rotates with 
a (nearly) Keplerian speed. The figure shows that the high-density Keplerian rotating disk 
(black line) is enclosed by a low-density pseudo disk. The pse udo-disk is formed by b oth the 
cloud rotation and magnetic field in the gas collapsing phase (IMachida et al. Il2005ai ). How- 
ever, the pseudo-disk gradually disappears in the main accretion phase, because its rotation 
speed does not reach the Keplerian velocity. Figures [5]a-c show that the driving point of the 
outflow radially expands with time, because the Keplerian rotating disk, which drives the 
outflow, extends in the radial direction with time. The outflow disappears in ~ 1.2 x 10 5 yr 
after the cloud collapse begins. In Figure EH, there is a sufficiently thin Keplerian rotation 
disk with a size of ~ 200 AU, while no outflow is driven by the circumstellar disk at this 
epoch. This is because by this epoch almost all the cloud mass already has fallen onto either 
the circumstellar disk or the protostar. The outflow is powered by the gas accreting energy 
onto the circumstellar disk. 

To investigate the evolution of the infalling envelope, the density distribution on the 
y = plane at the same epochs as in Figures [3] and [5] but for the initial cloud scale (~ 10 5 AU) 
is plotted in Figure |6j In each panel, the boundary between the gravitational sphere and 
static interstellar medium (r = 2tbe) is plotted by a white broken line; gas only inside 
the boundary can fall toward the center of the cloud. The figure shows that the cloud 
density at a larger radius gradually decreases with time, and the gas is concentrated onto 
the equatorial plane. The gas is also distributed along the z-axis in Figures |6f> and c. 
This mass distribution is caused by the mass ejection from the circumstellar disk by the 
protostellar outflow. Figure EH shows that gas in the infalling envelope is depleted and the 
mass ejection by the protostellar outflow has already stopped. It is natural that the elapsed 
time at Figure EH (t = 12 tff, c ) is longer than the freefall timescale (tff,b = 9.1 tff, c ) at the 
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cloud boundary, and thus the gas located near the cloud boundary at the initial state has 
already fallen onto the circumstellar disk by this epoch. It is expected that, at this epoch, 
magnetic braking becomes ineffective because the (less massive) infalling envelope cannot 
brake the circumstellar disk. Therefore, a large-size Keplerian rotating disk can appear at 
this epoch. 

Note that we may overestimate the mass of the infalling envelope, because we defined 
it as the mass of the whole cloud except for the Keplerian rotating disk. Thus, the infalling 
envelope includes the pseudo disk in our definition. The infalling envelope has a mass of 
M env = 0.23 M Q (i.e., 23% of the initial cloud mass) at t — 1.25 x 10 5 yr, as described in 
Figure |3]d In reality, however, a very small amount of the infalling gas remains in the region 
above the circumstellar disk as seen in Figure IHK 

The circumstellar disk is connected to the infalling envelope through the magnetic field 
lines. Thus, a slowly rotating massive infalling envelope can brake a rapidly rotating less 
massive disk through the magnetic field line. Thus, the magnetic field lines are strongly 
twisted when the infalling envelope is more massive than the circumstellar disk or when 
magnetic braking is effective. The difference between the rotation speeds of the infalling 
envelope and the circumstellar disk generates a strong toroidal field. Instead, when the 
infalling envelope is much less massive than the circumstellar disk, the infalling envelope 
cannot brake the circumstellar disk, and a strong toroidal field is not generated. In such a 
case, it is expected that the magnetic field has a nearly force free configuration. Figure [7] 
shows the ratio of the toroidal field to the poloidal field (B t /B p ) at the same epochs as in 
Figure [3j The toroidal field is greater than the poloidal field (toroidal dominated region; 
B r > Bp) inside the black line. Figure [7]a shows that the toroidal dominated region is 
distributed along the z-axis in the early main accretion phase, in which the torsional Alfven 
wave generated by the rotation of the circumstellar disk is propagated along the z-axis. 

In Figure [7]&, the toroidal field is much greater than the poloidal field in the region above 
the circumstellar disk (i.e., in the region of r c ^ 100 — 200 AU, where r c is the cylindrical 
radius). This indicates that magnetic braking is effective and the angular momentum of 
the circumstellar disk is transferred into the infalling envelope. At this epoch, strong mass 
ejection from the circumstellar disk also takes place as seen in Figure resulting from the 
generation of the strong toroidal field. The toroidal field weakens for t ^ 10tfj iC as seen in 
Figure [7]c. In Figure [7]d, the poloidal field dominates the toroidal field in the whole region 
of the infalling envelope. By this epoch, a large fraction of the mass in the initial cloud has 
already fallen onto either the circumstellar disk or the protostar. At this epoch, the mass 
of the infalling envelope is much less than the circumstellar disk mass (Fig. [6^). Thus, it 
is expected that, at this epoch, the magnetic field is rarely twisted and magnetic braking is 
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not so effective. 

Figure E] shows the radial distribution of the angular momentum for model 1 at the same 
epochs as in Figure E] against the radially integrated mass. For reference, the distribution 
of the angular momentum for unmagnetized model (model 5) at almost the same epoch as 
Figure [3]a is also plotted in this figure. At each radius, the angular velocity and mass are 
integrated in the gravitationally collapsing cloud (r < 2 r BE ) as 



where the protostellar mass M ps falling into the sink is added to estimate the mass. Note 
that, in Figure [8j the cloud mass does not reach M = 1M Q except for unmagnetized model 
(dotted line) because some fraction of the cloud mass (~ 10%) is ejected from the collapsing 
cloud (or from the region of r < 2r BE ) by the protostellar outflow. The figure indicates 
that the angular momentum of the whole cloud (M ~ 1M ) decreases gradually with time 
for model 1, while it is well conserved in the whole cloud for model 5. For model 1, the 
angular momentum of the cloud can be transferred into the interstellar medium by magnetic 
braking. However, the rate of reduction of the angular momentum is not very large during 
the calculation. The initial cloud has an angular momentum of J = 5.5 x 10 53 gcm 2 s _1 , 
while the collapsing cloud has J = 2.1 x 10 53 gcm 2 s _1 as the whole at the end of the 
calculation. Thus, about half of the angular momentum in the initial cloud is transferred 
into the interstellar medium. This is natural because the timescale of magnetic braking by the 
interstellar medium is longer than our calculation time (discussed in §1]). In addition, because 
the protostellar outflow driven by the circumstellar disk propagates into the interstellar 
medium with a certain amount of angular momentum, the protostellar outflow also transfers 
the angular momentum into the interstellar medium. 

Figure [H] also indicates that the angular momentum inside the disk (i.e., the region 
on the left side of the diamond symbol) is considerably smaller than that at the initial 
state. This difference of the angular momentum between the initial state and collapsing 
cloud clearly indicates that the angular momentum is effectively transferred outward in the 
collapsing cloud. In the early main accretion phase, the angular momentum is expected to 
be transferred by magnetic effects such as magnetic braking and protostellar outflow, since 
the circumstellar disk has an axisymmetric structure and remains clear axisymmetricity 
(Figs. [3]a and b), and hence, gravitational torque seems to be less important in magnetized 
models. In the later main accretion phase, since non axisymmetric patterns appear in the 
circumstellar disk (Figs. [3]c and d), the angular momentum is expected to be transferred 
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also by gravitational torque due to the non axisymmetric perturbation, in addition to the 
magnetic effects. 

3.2. Disk Formation in Clouds with Different Parameters 

In §3.1[ we investigated the disk formation for a typical model (model 1). However, 
it is expected that the formation process of the circumstellar disk and its size and mass 
depend on the initial cloud properties such as the magnetic field strength, rotation rate, and 
cloud mass (or size). In this subsection, we show the disk formation in clouds with different 
parameters. 

3.2.1. Different Magnetic Field Strengths 

To investigate the disk formation in clouds with initially different magnetic field strengths, 
we plotted the disk radii (a), disk and protostellar masses (b), disk-to-envelope mass ratios 
(c) and the masses of the infalling envelope (d) for models 1, 3, 4 and 5 against the elapsed 
time t c in Figure [91 The disk is defined as the region of R v < 1 (see, §3.1ft . Thus, in our 
definition, the disk (or the circumstellar disk) has a (nearly) Keplerian rotation. Note that 
the elapsed time t c , which is defined as the elapsed time after the circumstellar disk forma- 
tion, differs from t, which is defined as the elapsed time after the cloud collapse begins (or 
the calculation begins). They are related as 

t = t + t c , (18) 

where to is the time from the beginning of the calculation until the disk formation, and is 
listed in Table [2j In this subsection, we use t c in place of t to compare the disk evolution 
from the time just after the disk formation for different models. For models in Figure [91 
model 3 has a three times larger initial mass-to-flux ratio (A = 3) than model 1. Model 4 
has the same mass-to- flux ratio (A = 1) as model 1, but the cloud evolution is calculated 
with the ideal MHD approximation. The initially unmagnetized cloud (A = oo) is adopted 
for model 5. The elapsed time from the beginning of the calculation (t), disk radius {rf), 
mass of the circumstellar disk (Mdi S k) in each model at the end of the calculation are also 
listed in Table 

Figure [Hh indicates that the disk radius for an unmagnetized cloud (model 5) is much 
larger than the radii for magnetized clouds (models 1, 3, 4). In other words, in a magnetized 
cloud, the circumstellar disk is a relatively smaller that that in an unmagnetized cloud. 
This indicates that the magnetic field suppresses the disk formation, as already described 
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in past studies (Mellon fc Lill2008l . 120091 : Puffin fc Pudritzl 120091 : iHennebelle fc Ciardil 120091 ) . 
At t c ~ 10 5 yr, which almost corresponds to the freefall timescale at the cloud boundary, the 
circumstellar disk for model 5 reaches ~ 1200 AU, which corresponds to about 20% of the 
centrifugal radius of the whole cloud (see, Tabled]). For model 5, the increase and decrease of 
the disk radius for t > tg tC ~ 10 4 yr is attributed to the development of the spiral structure. 
As seen in Figure [9]fr, because the circumstellar disk mass is greater than the protostellar 
mass in the main accreti on phase, the disk is unstable against its gravity and frequently 



shows a spiral structure (IMachida et al.l l2010al ). The spiral structure effectively transfers 



the angular momentum outward and promotes a rapid mass accretion from the circumstellar 
disk onto the protostar. After the rapid mass accretion, the mass and size of circumstellar 
disk temporally decreases and the spiral structure weakens. For model 5, the circumstellar 
disk mass is greater than the mass of the infalling envelope for t c ^ 3 x 10 4 yr (Fig. [9]c). In 
addition, over 90% of the initial cloud mass already falls either onto the circumstellar disk or 
protostar at t c ~ 10 5 yr (Fig. [9]<f). Thus, it is expected that the main accretion phase almost 
ends by this epoch, and the circumstellar disk grows only a little by gas accretion after this 
epoch. 

On the other hand, at t c ~ 10 5 yr, the disk radii in magnetized clouds are about 100 — 
200 AU, which is considerably smaller than that of the unmagnetized cloud (~ 1200 AU). The 
large difference of disk radii between magnetized and unmagnetized models indicates that a 
large fraction of the angular momentum is transferred by magnetic effects. The difference 
of disk radii is caused by the different efficiencies of the angular momentum transfer among 
models. 

In Figure [9]a, resistive MHD models (models 1 and 3) show a similar trend for the 
evolution of the disk radius, in which model 3 has a three times larger mass-to-flux ratio 
than model 1 at the initial state. For these models, the disk radii gradually increases with 
time for t < (2 — 3) x 10 4 yr. Reflecting the difference in the initial mass-to-flux ratio, 
at the same epoch during the early main accretion phase [t < (2 — 3) x 10 4 yr], model 3 
which has an initially weaker magnetic field, has a slightly larger disk than model 1. Then, 
for t > (2 — 3) x 10 4 yr, the disk radii for these models increase with strong oscillation, 
which is caused by the spiral structure appearing in the circumstellar disk. Figure |H]fr shows 
that the circumstellar disk is more massive than the protostar by the end of the calculation 
even in magnetized clouds. In addition, the protostellar mass increases intermittently for 
t c ^ 2 x 10 4 yr. This is because the spiral structure, which appears in the circumstellar disk 
for tc ^ 10 4 yr (see, Fig. [3]c and d), unsteadily transfers the angular momentum outward 
and the mass accretion onto the protostar is temporarily amplified. Thus, during this phase 
(t ^ 10 4 yr), the angular momentum is transferred not only by the magnetic effects but also 
by the non axisymmetric perturbation in the circumstellar disk. 
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Finally, in resistive models (models 1 and 3), the disk radius rapidly increases for t c > 
(0.8 — 1) x 10 5 yr and reaches ~ 300 — 400 AU at the end of the calculation. Figure [9]c 
shows that the circumstellar disk mass is greater than the mass of the infalling envelope 
for t ^ tff 5 b — 10 5 yr. Note that since the infalling envelope includes the pseudo disk in 
our definition as described in §3.11 we somewhat overestimate the actual infalling envelope 
in Figure [9]c. After this epoch {t ^ 10 5 yr), it is expected that magnetic braking becomes 
ineffective, because the less massive infalling envelope can no longer brake the circumstellar 
disk. Thus, the circumstellar disk can rapidly grow without effective angular momentum 
transfer by the magnetic effects. As shown in Figure [9]g?, for models 1 and 3, over 80% of 
the initial cloud mass already falls either onto the circumstellar disk or the protostar by 
the end of the calculation. Note that a part of the initial cloud mass is ejected from the 
gravitationally collapsing cloud by the protostellar outflow as shown in Figure El At the end 
of the calculation, the circumstellar disk has a size of 290 AU for model 1 and 370 AU for 
model 3. Therefore, although the disk size in resistive MHD models is considerably smaller 
than that in the unmagnetized model, the Keplerian rotation disk with a size of > 100 AU 
can form even in strongly magnetized cloud. 



3.2.2. Resistive vs. Ideal MHD Models 

Figure [Hfc shows that the disk radius in the ideal MHD model (model 4) is somewhat 
smaller than that in the resistive MHD model (model 1) with the same parameters for the 
initial cloud. At end of the calculation (t ~ 2 x 10 5 yr), the circumstellar disk has a size of 
290 AU for model 1 and 120 AU for model 4. Magnetic braking is more effective in the ideal 
MHD model than in the resistive MHD model, because the magnetic field largely dissipates 
inside the circumstellar disk in the resistive MHD model. Thus, it is considered that the 
difference of the disk size between the ideal and resistive MHD models is caused by the 
difference in the strengths of the magnetic field in the circumstellar disk. 

Figure [TU] shows the plasma beta for models 1 (upper panel) and 4 (lower panel) on the 
y = plane at almost the same epoch after the calculation begins. The contour of (3 P = 1 
is shown by a white thick line, outside of which the magnetic energy is greater than the 
thermal energy. The circumstellar disk exists in the region of (3 P ^ 1 (i.e., inside the white 
line). Because the circumstellar disk has a high density (or large thermal pressure), the 
thermal energy is greater than the magnetic energy (i.e., (3 P > 1) inside the circumstellar 
disk. Figure [10] shows that, for model 1, the plasma beta is extremely high as /3 P ^ 10 5 in 
the region of \r\ <^ 50 AU inside the circumstellar disk, because this region has a density of 
n ~ 10 12 cm -3 (see Fig. [5]c) and the Ohmic dissipation is effective. Figure [10] lower panel 
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shows that the circumstellar disk for the ideal MHD model is smaller than that for the 
resistive MHD model. In the circumstellar disk, the plasma beta for the ideal MHD model 
is about two or three orders of magnitude smaller than that for the resistive MHD model. 
Despite the large difference in the magnetic energy in the circumstellar disk, the difference 
in the final size of the circumstellar disk is not very large: the disk has a size of 290 AU for 
model 1, and 120 AU for model 4. In addition, the evolution of the masses of circumstellar 
disk, protostar, and infalling envelope for the ideal MHD model is almost the same as that 
for the resistive MHD model, as shown in Figure [9]& and c. Thus, even in the ideal MHD 
model, magnetic braking becomes ineffective as the mass of the infalling envelope is depleted, 
and a circumstellar disk with a size of ~ 100 AU can form. 



3.2.3. Effect of Initial Cloud Rotation 

To investigate the effect of the initial cloud rotation on the circumstellar disk formation, 
the disk radii and circumstellar disk and protostellar masses for models 1, 6 and 7 are plotted 
in Figure [TT1 Model 6 has three times the initial angular velocity as model 1, while model 7 
has one third the initial angular velocity as model 1. 

Figure [TT] upper panel shows that at the same epoch after the circumstellar disk forma- 
tion, the model with rapid initial rotation has a relatively large size disk. In addition, the 
disk size begins to oscillate from the earlier evolution phase in the cloud with faster rotation. 
For example, model 6 begins to show an oscillation at t c ~ 2 x 10 4 yr, while model 7 shows 
it for t c ^ 5 x 10 4 yr. The size oscillation occurs after the spiral structure develops in the 
circumstellar dis k. The spiral structure appears after the circumstellar disk grows sufficiently 



in size and mass (IMachida et al.ll2010al ). Figure [T2l shows the density distribution for models 
6 and 7 on the equatorial plane after the circumstellar disk becomes more massive than 
the infalling envelope. The figure indicates that the circumstellar disks show a clear spiral 
structure with a size of 300 — 400 AU. 

Figure [TT] lower panel shows that the cloud with faster initial rotation creates a more 
massive circumstellar disk. Thus, in such a cloud, the spiral structure and size oscillation 
appear in an earlier phase of the disk evolution. In addition, the cloud with faster initial 
rotation has a less massive protostar before the size oscillation because mass accretion from 
the circumstellar disk onto the protostar is suppressed by a strong centrifugal force in such 
a model. After the size oscillation occurs, the protostar unsteadily acquires its mass from 
the circumstellar disk, as described in §3.11 Finally, the disk radius in each model increases 
rapidly for t c £ (0.7 - 0.9) x 10 4 yr to reach ~ 300 - 400 AU (see, Table W)- This is because, 
by this epoch, the circumstellar disk mass becomes greater than the protostellar mass, and 



- 19 - 



magnetic braking becomes ineffective. 

3.2.4- Effect of Initial Cloud Masses 

Figure [13] plots the disk radius and mass of the infalling envelope normalized by the 
initial cloud mass for models 1, 8, and 9 against the elapsed time after the circumstellar disk 
formation. Each of these models has a different initial cloud mass but the same A and 
The figure shows that the initially less massive cloud has a larger-size disk at the same epoch 
after the circumstellar disk formation. The disk radii for models 1 and 8 rapidly increase 
after the mass of the infalling envelope deceases to < 20 — 40% of the initial cloud mass. 
For model 8, which has an initial cloud mass of M c \ = 0.3 M , the mass of the infalling 
envelope normalized by the initial cloud mass decreases to ~ 0.2 at t c ~ 10 4 yr, and the 
circumstellar disk radius rapidly increases for t c ^ 10 4 yr and reaches ~ 220 AU by the 
end of the calculation. On the other hand, for model 9, which has an initial cloud mass 
of M c i = 10 Mq, the size of circumstellar disk does not exceed 30 AU by the end of the 
calculation (t c ~ 10 5 yr). Thus, the disk size for model 9 is less than one-tenth of that for 
model 1 at t c ~ 10 5 yr. Note that, for model 9, because the freefall timescale at the cloud 
boundary is as large as ~ 10 6 yr, we cannot calculate the cloud evolution until the infalling 
envelope is depleted. 

It is expected that the difference in the disk size for models with different initial cloud 
mass is caused by the remaining mass of the infalling envelope. As described in ^TJ the 
massive infalling envelope brakes the circumstellar disk through the magnetic field lines. 
The initial cloud mass and freefall timescale at the center of the cloud or at the cloud 
boundary are different in models 1, 8, and 9. Thus, at the same epoch, a model with an 
initially massive cloud has a massive infalling envelope that can brake the circumstellar disk 
and suppresses the growth of the circumstellar disk by magnetic braking. It is difficult to 
quantitatively relate the disk size and the mass of the infalling envelope among these models, 
because the initial distribution of the mass-to-flux ratio (see Fig. [2]) and angular velocities 
(see Table [1]) are different in these models. However, Figure [1] indicates that after the mass 
of the infalling envelope decreases sufficiently, the effect of magnetic braking weakens and 
a large-size circumstellar disk can form even in a strongly magnetized cloud. Conversely, 
large-disk formation is suppressed by magnetic braking when a massive infalling envelope 
remains around the protostar. 
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Discussion and Conclusion 



In this paper, to investigate the formation of the circumstellar disk in a realistic setting, 
we set a dense molecular cloud core embedded in a low-density interstellar medium in a 
large simulation box with size 16 times the radius of the molecular cloud core, in which 
a uniform magnetic field parallel to the rotation axis is adopted. In the gravitationally 
collapsing cloud, it is expected that the disk formation is suppressed by magnetic braking 
which effectivel y transfers the angular momentum from the circumstellar disk, as describe d 
in past studies jMellon & Lill2008l . bood : Irlennebelle & Ciardil[2009[ buffin & Pudritzll2009} ). 
In reality, our results indicate that the circumstellar disk formed in a magnetized cloud is 
much smaller than that in an unmagnetized cloud (see §3.2. ip . However, we also showed 
that the circumstellar disk exceeding ~ 100 AU can form even in a strongly magnetized 
cloud (A ~ 1), which contradicts previous studies. The past studies predicted that magnetic 
braking suppresses the disk formation and limits the size of the circumstellar disk within 
~ 10 AU. This inconsistency is attributed to the initial setting of the host cloud and the 
calculation time after the circumstellar disk formation. In this study, we calculated the 
cloud evolution until almost all gas inside the initial cloud falls onto either the circumstellar 
disk or the protostar. In other words, we execute the calculation until the freefall timescale at 
the cloud boundary is exceeded. As described in §[TJ as the cloud mass is depleted, magnetic 
braking is expected to be ineffective because the less massive infalling envelope cannot brake 
the circumstellar disk. In §31 we numerically showed that magnetic braking weakens with 
time and a large-size circumstellar disk can form. In this section, we analytically discuss the 
timescale of magnetic braking and disk formation. 

The angular momentum of the cloud or disk is transferred into the external medium 
through a torsional Alfven wave. Magnetic braking becomes effective when the angular 
momentum transferred into the external medium is comparable to that of the rotator (e.g., 
the molecular cloud or disk). In other words, magnetic braking becomes effective when the 
moment of inertia of the gas through which the Alfven waves have propagated becomes 
comparable to the moment of inertia of the cloud or disk. With this condition, the magnetic 
braking timescale can be described as 
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where z d , v^ >eK , pd and p ex are the half thickness of the cloud (or the disk), Alfven veloc- 
ity in the ex ternal medium, cloud ( or disk) density, and densi t y of the external medium, 
respe ctively (jMestel fc Spitzerl Il956l : iNakano &: Tademarul 1 19721 ; iMouschovias &: Paleologou 
1980j). Thus, magnetic braking becomes effective when a typical timescale of the system (r) 
is longer than the magnetic braking timescale r > rn. With equation (fT9l) . the condition for 
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magnetic braking becoming effective is described as 

rv A , cx p cx ^ z d p d . (20) 

The left-hand side of equation (120!) corresponds to the line density of the external medium 
in the range at which the Alfven wave can reach in r, while the right-hand side corresponds 
to the line density of the cloud or disk. Equation ( ]20]) indicates that magnetic braking 
is effective only when the mass of the external medium is greater than the mass of the 
rotator. In other words, magnetic braking does not take place when sufficient mass does 
not exit around the rotator. This is natural because a less massive external medium cannot 
brake a massive object. This condition can be also understood from equation f JT9|) . When 
an uniform magnetic field is assumed, equation f[T^]) indicates that the magnetic braking 
timescale is proportional to t\\ oc p cxt . Note that the Alfven speed is proportional to p ext ' . 
Thus, the magnetic braking timescale becomes longer as the external density decreases. For 
example, when the external density is extremely low, magnetic braking does not take place. 

We estimate the magnetic braking timescale for model 1 (initially the strongest magnetic 
field model). First, we estimate the magnetic-braking timescale of the molecular cloud core 
embedded in the interstellar medium. For model 1, the density contrast between the center 
of the cloud core and the external medium is 82 at the initial state, and Alfven speed in 
the external mediums is t>A,ext = 1.3kms _1 . Thus, with equation (TT9]) . the magnetic braking 
timescale becomes th ~ 8 x 10 5 yr, which is much longer than the freefall timescale of the 
cloud (tff )C = 1.04 x 10 4 yr) or our calculation time (~ 10 5 yr). Thus, it is expected that the 
cloud rarely loses its angular momentum by magnetic braking of the interstellar medium. 
In §3.11 we showed that the total angular momentum of the cloud core is rarely decreased 
through the calculation. 

Next, we estimate the magnetic-braking timescale of the circumstellar disk. First, we 
assume that the circumstellar disk exists only in the low-density interstellar medium (thus, 
not including the infalling envelope) to simply estimate the effect of the low-density interstel- 
lar medium outside the gravitational sphere. The density contrast between the circumstellar 
disk (n ~ 10 10 cm -3 ) and interstellar medium (3 x 10 4 cm~ 3 ) is about 3 x 10 6 , and the cir- 
cumstellar disk has a scale height of ~ 10 AU. With these parameters, the magnetic braking 
timescale of the circumstellar disk is tii ~ 1.1 X 10 8 yr. Thus, the interstellar medium cannot 
brake the circumstellar disk in our setting. Note that, in this study, we adopted a somewhat 
higher density of the interstellar medium (3 x 10 4 cm -3 ) than the observed value of the inter- 
stellar medium to limit the Alfven speed in the interstellar medium, because it is difficult to 
calculate the cloud evolution over a long time with a large Alfven speed. Thus, the magnetic 
braking timescale in the actual interstellar medium is longer than in our setting with a lower 

— 1/2 

external density, because it is proportional to p ext . In summary, magnetic braking from the 
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low-density interstellar medium is completely ignored in the formation and evolution of the 
circumstellar disk. 

On the other hand, the infalling envelope in the region of r < 2r BE can brake the 
circumstellar disk in the main accretion phase. The torsional Alfven wave generated by the 
center of the cloud can reach the cloud boundary (r = 5400 AU) in ~ 2 x 10 5 yr (~ 2tff jb ), 
where the Alfven speed of the external medium VA,ext = 1.3 km s" 1 is adopted. Note that, 
in reality, because the magnetic field is amplified, the torsional Alfven wave can reach the 
cloud boundary in several freefall timescales. Thus, when the infalling envelope is more 
massive than the circumstellar disk, the angular momentum of the circumstellar disk is 
effectively transferred by magnetic braking. In the main accretion phase, however, the mass 
of the infalling envelope decreases with time. Thus, it is expected that magnetic braking 
gradually becomes ineffective with time. However, it is difficult to estimate the magnetic 
braking timescale of the infalling envelope quantitatively, because the infalling envelope 
drastically changes its density distribution in a several free-fall timescale that corresponds 
to the timescale for the circumstellar disk formation. 

Thus, the disk formation process and final size of the circumstellar disk cannot be inves- 
tigated using a simple analytical framework. Only a less-massive and small-size circumstellar 
disk may appear in a magnetized cloud when magnetic braking significantly transfers the 
angular momentum by the end of the main accretion phase. In this study, however, we 
calculated the cloud evolution by the end of the main accretion phase and showed that the 
circumstellar disk exceeding ~ 100 AU can form even in a strongly magnetized cloud. In a 
magnetized cloud, magnetic braking becomes ineffective and the disk size increases rapidly 
in the later main accretion phase. It is natural that magnetic braking only redistributes the 
angular momentum from the circumstellar disk into the infalling envelope. Because the in- 
falling envelope finally falls onto the circumstellar disk, the excess angular momentum gives 
rise to a rotation supporting the disk in the later main accretion phase. Instead of magnetic 
braking, only the protostellar outflow can transfer the angular momentum to the interstellar 
medium, because the gas with a certain amount of angular momentum is directly ejected 
into the interstellar medium from the circumstellar disk by the protostellar outflow. Because 
the protostellar outflow is powered by mass accretion onto the circumstellar disk, the out- 
flow disappears in the later main accretion phase, which also corresponds to rapid growth 
of the circumstellar disk. Thus, the circumstellar disk begins to grow after the protostellar 
outflow sufficiently weakens or disappears. In addition, in the later main accretion phase, 
the spiral structure appearing in the circumstellar disk can contribute to the redistribution 
of the angular momentum in the circumstellar disk. Thus, it is difficult to specify which 
mechanism is effective for the angular momentum transfer in a magnetized collapsing cloud. 
However, the large-size disk comparable to observations can form in a strongly magnetized 
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cloud. 

In this study, we mainly investigated the evolution of the collapsing cloud with a mass of 
1 Mq in the initial state and found that the circumstellar disk with a size of several hundred 
AU appears in the later main accretion phase. On the other hand, we cannot calculate 
the cloud evolution for an initially massive cloud until the infalling envelope disappears. 
However, a massive cloud has a larger angular momentum (or larger centrifugal radius). 
Thus, a circumstellar disk much larger than ~ 100 AU may form in an initially massive 
cloud. 

Finally, we comment on the geometry of the initial magnetic field. In this study, we 
assumed the magnetic field lin e s to b e parallel to the rotation axis. With ideal MHD cal- 



culation, iHennebelle fc Ciardil (j2009f ) already reported the formation of the circumstellar 
disk when the magnetic field lines are inclined from the rotation axis. Thus, we need to 
investigate the inclined field case with non-ideal MHD effects in the future. 

Numerical computations were carried out on NEC SX-9 at Center for Computational 
Astrophysics, CfCA, of National Astronomical Observatory of Japan, and NEC SX-8 at the 
Yukawa Institute Computer Facility. This work was supported by the Grants-in-Aid from 
MEXT (20540238, 21740136). 
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Table 1: Model parameters 



Model 


n Cj0 [cm- 3 


r d [AU] 


Mci [M ] 


a 


UJ 


A 


B [mG] 


o [s- 1 ] 


Po 


Ra [au] R/r 


1 


3 x 10 6 


5400 


1 


1.0 


0.1 


1 


176 


2.3 x 10" 13 


0.03 


2150 


R 


2 


3 x 10 6 


5400 


1 


0.3 


0.1 


2 


92 


2.3 x 10~ 13 


0.03 


2150 


R 


3 


3 x 10 6 


5400 


1 


0.1 


0.1 


3 


58 


2.3 x 10~ 13 


0.03 


2150 


R 


4 


3 x 10 6 


5400 


1 


1.0 


0.1 


1 


176 


2.3 x 10~ 13 


0.03 


2150 


I 


5 


3 x 10 6 


5400 


1 





0.1 


oo 





2.3 x 10~ 13 


0.03 


2150 


R 


6 


3 x 10 6 


5400 


1 


1.0 


0.15 


1 


176 


4.7 x 10~ 13 


0.04 


4800 


R 


7 


3 x 10 6 


5400 


1 


1.0 


0.05 


1 


176 


1.2 x 10~ 13 


0.007 


540 


R 


8 


3 x 10 7 


1740 


0.3 


1.0 


0.1 


1 


566 


7.4 x 10~ 13 


0.03 


680 


R 


9 


3 x 10 4 


54000 


10 


1.0 


0.1 


1 


17.6 


2.3 x 10~ 14 


0.03 


21500 


R 



* R: Resistive MHD model, I: Ideal MHD model 



Table 2: Calculation results 



Model 


t [tff.c] 


to [*ff,c] 


r d [AU] 


M disk [M ] 


1 


12.8 


3.39 


290 


0.49 


2 


12.3 


3.09 


440 


0.39 


3 


13.2 


2.95 


370 


0.42 


4 


11.2 


3.40 


120 


0.31 


5 


11.9 


2.87 


1200 


0.91 


6 


11.3 


3.57 


390 


0.50 


7 


11.4 


3.38 


410 


0.52 


8 


14.2 


3.59 


220 


0.12 


9 


4.41 


3.37 


29 


0.05 
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r [AU] 

Fig. 1. — Density distribution (left axis) and enclosed mass (right axis) are plotted against 
the initial cloud radius for models with a central density n C) o = 3 x 10 6 cm -3 (thick line) and 
3 x 10 7 cm -3 (thin line). Only the enclosed mass in the range of r < 2tbe is plotted. 
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Fig. 2. — Initial mass-to- flux ratios normalized by the critical value (A) for models with 
different initial magnetic field strengths a or different cloud masses are plotted against the 
initial cloud mass (left panel) and cloud radius (right panel) in the upper panels. The 
centrifugal radii (R c f) for different angular velocities oj or different cloud masses are plotted 
in the lower panels. 




Fig. 3. — Density distribution {color and contours) and velocity vectors (arrows) on the 
equatorial plane at four different epochs for model 1. The elapsed time t, protostellar mass 
M ps , circumstellar disk mass M disk , and envelope mass M env are listed in each panel. 
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Fig. 4. — (a) Azimuthal velocity normalized by the Kepler velocity, (b) ratio of the radial 
velocity to the azimuthal velocity, (c) surface density and (d) plasma beta /3 p for model 1 at 
the same epochs as Figs. [3] are plotted against the cloud radius. Each quantity is azimuthally 
averaged. 
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(a) t =3.758 tfi.o = 3.94x1 O a [yr] (b) t =7,307 t tf .o = 7.67x1 0* [y r] 
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Fig. 5. — Same as Figure |3]but on the y = plane. In each panel, the black line corresponds 
to the circumstellar disk. The boundary between outflowing (v r > 0) and inflowing (v r < 0) 
gas is described by the white line; inside this boundary, the gas is outflowing from the 
circumstellar disk. 
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Fig. 6. — Density distribution (color and contours) for model 1 in the region of z > on 
the y = plane at the same epochs as Figs. [3] and [5] except for panel a (initial state). The 
white broken line is the boundary between the gravitational sphere and interstellar medium; 
gravity (gas self-gravity and the gravity of the protostar) is imposed inside this boundary. 
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Fig. 7. — Ratio of the toroidal field to the poloidal field {color and contours) and magnetic 
field (arrows) for model 1 in the region of z > on the y = plane at the same epochs as 
Figures EJ [5] and |6j The toroidal component (B t ) is greater than the poloidal component 
(B p ) inside the black line. The white line corresponds to the contour of B t = 0.1B P . The 
arrows indicate the magnetic field direction on the y = plane. 
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Fig. 8. — Distribution of the angular momentum for model 1 at the same epochs as in 
Figure [3] against the cloud mass. The dashed line is the distribution of the angular momentum 
for model 5 at t = 1.67 x 10 4 yr. The initial distribution of the angular momentum is plotted 
by a thin solid line. The angular momentum and mass are integrated in the region of < r. 
The sum of the protostellar mass M ps and circumstellar disk mass of each epoch is 

plotted by the diamond symbol (o) on each line. The region on the left-hand side of the 
diamond symbol corresponds to the circumstellar disk region. 
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Fig. 9. — Disk radii (a), disk-to-envelope mass ratios (b), circumstellar disk and protostellar 
masses (c) and masses of the infalling envelope (d) for models 1, 3, 4 and 5 against the 
elapsed time after the circumstellar disk formation. The freefall timescale at the center of 
the cloud tg iC and at the cloud boundary tg t b are plotted in each panel. 
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Fig. 10. — Plasma beta (color) on the y = plane for the resistive (model 1; upper panel) 
and ideal (model 4; lower panel) MHD models. The thick white line is the contour of (3 P = 1. 
The arrows indicate the magnetic field direction on the y = plane. 
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Fig. 11.— Disk radii (upper panel) and disk and protostellar masses (lower panel) for 
models 1, 6 and 7 against the elapsed time after the circumstellar disk formation. The 
freefall timescale at the center of the cloud tg tC and at the cloud boundary tg^ are plotted 
in each panel. 
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Fig. 12. — Density distribution {color and contours) and velocity vectors (arrows) on the 
equatorial plane for models 6 and 7. The elapsed time t is described in each panel. 
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Fig. 13. — Disk radius (left axis) and mass of the infalling envelope normalized by the 
initial cloud mass are plotted against the elapsed time after the circumstellar disk formation 
for models 1, 8 and 9. The freefall timescale for each model are noted by the circle on the 
horizontal axis. 



